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Abstract 



I We have constructed, numerically, both regular and black hole static 

■ solutions to the simplest possible gravitating Yang-Mills-Higgs (YMH) in 
. 4p spacetime dimensions. The YMH systems consist of 2p— th power cur- 
vature fields without a Higgs potential. The gravitational systems consist 

■ of the 'Ricci scalar' of the p— th power of the Riemann curvature. In 4 

■ spacetime dimensions this is the usual Einstein-YMH (EYMH) studied in 
^ ' [1) 2], whose qualitative results we emulate exactly. 



1 Introduction 

Gravitating monopoles were studied intensively in [1, 2, 3], subsequent to the 
discovery of the finite energy regular solutions [4] to the gravitating Yang-Mills 
(EYM) system. These regular EYM solutions were soon extended to black 
holes [5]. EYM solutions in 3 -|- 1 dimensional spacetime are reviewed in [6] 
and an exhaustive analysis is given in [7]. Along with the introduction of the 
(dimensionful) Higgs field that leads to the construction of regular and black 
monopoles [1, 2, 3] in the gravitating EYM-Higgs (EYMH) system, the (dimen- 
sionful) cosmological constant was introduced. The case of positive cosmological 
constant is reviewed quite adequately in [6] and that of negative is in the review 
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[8]. The latter led to very interesting new features of the solutions and is, in some 
sense, an alternative to the Higgs field. 

EYM solutions in (higher) dimensions D + 1 {D > A) were considered only 
relatively recently. Here there arc two possibilities of identifying the D spacelike 
dimensions: Either as asymptotically Minkowski with S^~^ boundary, or, with 
boundary S^~'^'~^ x in which case the N codimensions are frozen as in the 
case of the 2;— coordinate of the Abrikosov-Nielsen-Olesen vortex. In the case 
of boundary for D = A it was found in [9] that the energy of the regular 
solution is infinite, and in [10] (with negative cosmological constant) it turned 
out that the energy of the black hole is also infinite. In addition, in [9] solutions 
with boundary S'^ x D — A with finite energy per unit length along R^ were 
constructed, the total energy remaining infinite. This is not surprising because 
the usual EYM system in D > A does not have the requisite scaling properties for 
there to exist finite energy solutions. This obstacle is circumvented by extending 
the definition of the EYM system to feature higher order YM curvature terms 
with the appropriate scaling. 

Suitably defined EYM systems featuring higher order YM curvature^ terms in 
aAl d = D + 1 spacctimc dimensions were studied in [11, 12] withoTit cosmological 
constant, and in [13, 14] with cosmological constant. These latter differ from their 
d— 3+1 dimensional EYM analogues [4, 6], which are unstable (sphaleronic) and 
have no gravity decoupling limits, in that some of them do have gravity decoupling 
limits which are stabilised by Chern-Pontryagin (instantonic) charges. In this 
respect these are more akin to the gravitating monopoles, and like them these 
systems exhibit one or more dimensionful constants that cannot be scaled away 
and hence parametrise the ensuing solutions. As a result they feature certain 
bifurcation properties like the monopoles but some of these exhibit in addition 
what were named conical singularities in [15], where this analysis was carried 
out. To date no gravitating monopoles in dimensions higher than d = 3 + 1 are 
constructed and it is the aim of the present work to do that, in the simplest 
class of YMH models defined in Ap spacetime dimensions. Before going into the 
details of this choice, we point out that this is a very interesting restriction since 
in our previous study of higher dimensional EYM solutions in [15], it was found 
that the qualitative features of these repeated modulo every 4— dimensions. Like 
in all the work on gravitating gauge fields quoted above, our solutions take into 
account the backreaction of gravity on the YMH fields. 

Before considering the possible choices of suitable higher dimensional YMH 
systems, it is in order to consider the particular interest in gravitating mono- 
poles in higher dimensions, a) on a technical level, and b) from the viewpoint of 

^Higher order Riemann curvature terms can also be introduced, but since all studies are 
in practice carried out for systems subject to symmetries, such terms of high enough order 
Riemann curvature vanish due to the symmetry imposed. Their inclusion is thus unnecessary. 
Any nonvanishing such terms included do not play an essential role, but rather result only in 
a quantitative difference. 
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applications. 

a) On a technical level, it is the only way which enables the construction 
of static solutions of EYM systems in spacetime dimensions higher than four, 
which describe a nonvanishing 'non Abelian electric' connection Aq. In this re- 
spect EYM solutions in the presence of a negative cosmological constant differ 
in spacetime dimension d — A, from those in d > 5. In d — A the presence of 
the negative cosmological constant results in such asymptotic properties, which 
in contrast to the asymptotically fiat case, allow for nonvanishing electric field. 
This is a consequence of the fact that the 'magnetic' connection of these solu- 
tions can [8], among other possibilities, behave as a half pure gauge iust like a 
monopole. We have verified that the corresponding EYM systems (with nega- 
tive cosmological constant) in > 5 support 'magnetic' connections that behave 
asymptotically as pure gauge, like their asymptotically fiat counterparts. Hence 
they do not support a nonvanishing electric field. 

When a Higgs field is introduced however (irrespective of the presence of 
a cosmological constant), the situation changes. Already in fiat space, EYMH 
systems automatically support dyonic solutions, e.g., in all even dimensional 
spacetimes for models employed in [16]. Clearly, when the gravitational force is 
switched on, the 'electric' field of the resulting EYMH solution will persist. 

b) Prom the viewpoint of applications, like their higher dimensional EYM 
analogues, they are expected to be relevant to various aspects of the study of 

branes. Some examples in the literature concern monopoles in string the- 
ory [17, 18], and self gravitating supersymmetric solitons in [19]. In the former 
examples [17, 18] only the Yang-Mills system in the absence of Higgs fields ap- 
pears, so that their higher dimensional extensions will feature the gauge field 
configurations described by the Yang-Mills hierarchy.^ The work of [19] on the 
other hand concerns monopoles of the usual YMH system in the spacelike sub- 
space. Thus its higher dimensional extensions will involve field configurations of 
the YM- Higgs hierarchies (see below), in other words higher dimensional mono- 
poles. What is more, is that the monopoles in [19] satisfy the selfduahty constraint 
in the spacelike dimensions, and that the higher dimensional extensions of the 
YMH systems employed here arc chosen precisely such that the corresponding 
extended Bogomol'nyi constraints are likewise in force. 

When choosing a YMH system in arbitrary (Euclidean) dimension D support- 
ing a monopole solution we are faced with a plethora models. The most efficient 
way of constructing these YMH models is via dimensional descent [21] from a YM 
system [20] on D + N (even) dimension. Integrating out the coordinates on the 
(compact) codimensions results in the dimensional residual YMH theory 
supporting monopoles, whose topological (monopole) charge is the descendent of 

^In [17, 18] instanton fields of the usual YM system are exploited, while their higher dimen- 
sional extensions would employ the corresponding instantons of the YM hierarchy introduced 
in [20]. 
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the |(-D + N)—th Chern-Pontryagin charge. We shall not dwell on the detailed 
properties of these D— dimensional monopolcs here, save to emphasise their most 
relevant feature pertinent to the present work. Because of the high degree of 
nonlinearity of the extended selfuality equations [20] in D + N dimensions, the 
descendent Bogomol'nyi equations in D dimensions are overdetermined [22] and 
in general cannot be saturated. It turns out that these Bogomol'nyi equations can 
be saturated only when the descent is over = 1 and N = D — 2 codimensions, 
the latter case being irrelevant here. Since the Euclidean dimensions D + N are 
only even, then the residual Euclidean space is D — 1 dimensional, restricting us 
to models in even spacetime dimensions d — {D — 1) + 1 only. 

In the present work we have restricted our attention to d — Ap dimensional 
spacetimes because the YM hierarchy in D + 1 = 4p dimensions is scale invariant 
and does not feature an additional dimensionful constant.^ This results in the 
simplest possible residual YMH models, keeping a tight analogy with the 't Hooft- 
Polyakov monopole (in the BPS hmit). Like the latter, these models feature 
only one dimensional constant, namely the Higgs VEV with inverse dimension 
of a length. A more direct construction, with the sole criterion of achieving a 
topological lower bound and not applying dimensional descent, was employed in 
[24, 25, 16], but these models always feature an extra dimensional constant in 
addition to the Higgs VEV, considerably complicating the numerical analysis in 
the gravitating case. Thus we have eschewed such models here. 

Finally, we state the action densities to be employed. It is convenient to 
express the connection and the Higgs field in terms of the chiral so(4p) represen- 
tation matrices 

'^nu ^ {'^ij,'^i,D+i) , // = 0,i, i^l,2,...,D , (1) 

such that the connection takes its values in the algebra of S0{4:p— 1), namely 
Ejj, and the Higgs field is 

$ = 0^E,,,,+i. (2) 

In terms of the curvature F = F{2) of A and the covariant derivative D^, the 
action densities of this hierarchy of YMH models in 4p spacetime are defined as 



C(4P) _ rr. 

'-'matter — J-l 



F{2pf r7{F{2p - 2) A D^f 



2{2p)\ ' 2(2p-l)! 



(3) 



the 2A;— form F{2k) being the A;— fold totally antisymmetrised product of F = 
F{2). The flat space static dyons of these systems were constructed in [16]. It 
is the monopole solutions in [16] that are gravitated in the present work, in the 
spirit of [1, 2]. 



•^YM systems on R^'^^^ supporting instantons involve at least one additional dimensionful 
constant (see e.g., [23]) whose dimensional descendent therefore will feature one more dimen- 
sional constant in addition to the Higgs VEV. 
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A pertinent comment at this point is that had we chosen the Higgs kinetic term 
in (3) to be {F{2k — 2) A D^y, with k ^ p in 2{p + k) dimensions, the sohitions 
would still be topologically stable [16] but now a new dimensional constant will 
appear in (3). In this respect such systems would be more akin to the higher 
dimensional models studied in [15] and their gravitating solutions could then 
exhibit conical fixed points. The fixed point analysis for such selfgravitating 
YMH has not been carried out to date and it promises to be appreciably more 
involved than the corresponding analysis in [15], basically because now there is 
a Higgs field function in addition to the gauge field function. We have eschewed 
this choice here keeping strictly the analogy with the gravitating monopoles in 
3 + 1 dimensions, and encounter no solutions featuring conical fixed points. We 
will return to this question elsewhere. 

To complete the definition of the gravitating YMH systems, the gravitational 
Lagrangian must be specified. Restricting to Levi-Civita connections, the hier- 
archy of Einstein systems in d = 2{p + q) dimenional spacetimes are 



where E{2p) = *e(2g)* is the double-Hodge dual of e(2g), the 2^— fold anti- 
symmetrised product of the Vielbein fields, e in (4) is the determinant of the 
Vielbein and R{2p) is the 2p— form Riemann curvature. The p = 1 member of 
the hierarchy (4) is the Einstein-Hilbert Lagrangian in d — 2{q + 1), p — 2 the 
Gauss-Bonnet system in d = 2{q + 2), etc. Our choice of Lagrangian employed 
for gravitating the YMH system (3) for given p is the p = q = ^ member of 
(4), like in the p = 1 case studied in [1, 2]. One can of course choose the p = 1 
member of (4) with all YMH actions (3) but in that case the relative dimension- 
ality of the YMH and gravitational Lagrangians would be [L^^^^"^)] instead of 
[L^*'], which then would match with the p — q — 1 case studied in [1, 2], only for 
p = 1. This would result in the dilution of the analogy with the latter at least 
on the quantitative level. This choice was made previously in the study of the 
EYM systems in 4p dimensional spacetimes in [26], where it resulted in the tight 
similarity in the qualitative properties of the solutions of the p— th YM system 
gravitated with the p{— q)—th Einstein systems in all 4p dimensions. 

The hierarchy of monopoles to be studied here is that of static and spherically 
symmetric (in Ap — 1 space dimensions) solutions to the equations arising from 
the action densities 



2 Ansatz, Action, and Differential Equations 

Using the spherically symmetric metric Ansatz in 'Schwarzschild' coordinates 



e%,,) = (E(2p),i?(2p)), 



(4) 




(5) 



ds^ = A\T)ii{r)de 



dr^ 



(6) 



/i(r) 



{d-2) > 
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together with a generahsed 't Hooft-Polyakov Ansatz for the magnetic monopole 
^ = ^,,4(1-^0), cl>^^-h{r). (7) 



we obtain the reduced action 



j[\/Jiy (~^p^o + f,SM), (8) 
with Hp = Kp/22(p-i), fp = Tp/(2p)!, and 

Sa^-ydrA^^{r'-^^-\l-t,r), (9a) 



2p 



(9b) 



where W = {w^ - if and Hp = (w^ - lyP'^h. 

A Derrick-type scahng argument shows that static finite energy solutions of 
the field eqs. derived form the action (8) can only exist in spacetime dimensions 
d with 

2p+ 1 < d < 4p+ 1 , (10) 

and indeed in [26] the EYM system has been studied for all values in this range. 
The presence of the Higgs field, however, requires d = 4p. 

The action (8) depends on three dimensionful parameters, kp, fp, and the vac- 
uum expectation value 77 = lim^^oo h{r) of the Higgs field, with the dimensionless 
ratio 

Since an overall factor in front of the action has no effect on the eqs. of motion, 
we can rescale r such that either 

S = const. • {Sg + oi^^Sm) , = 1 , (12) 

resulting in a-dependent eqs. of motion with a-independent boundary conditions, 
or 

5" = const. • {Sg + Sm) , ^ = a , (13) 

resulting in a-independent eqs. of motion with an a-dependent boundary condi- 
tion. 

These two formulations are clearly equivalent as long as a 7^ 0. In the 
limit a ^ they are, however, quite different. The first alternative describes 
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a YMH system in a fixed gravitational background, either flat (no gravity) or a 
Schwarzschild type black hole. The second alternative describes the (self gravi- 
tating) EYM system of [26]. 

In the following we will mostly use the second formulation, but will occasion- 
ally refer to the first formulation as 'unsealed variables'. 

The dependence of the action (8) on A{r) suggests to define the 'mass function' 

M(r) = ^^±i^/-^--^(l-/.)^ (14) 

that has a finite limit M = liuir^^ M{r) for asymptotically fiat solutions with 
finite energy. 

For the Reissner- Nordstrom solution with w = and h = a the mass function 
satisfies 

^ = ^-y-\ 4p + 1 - d)r-(^-+-'^) , (15a) 
ar 2p 

M{r) = M- l:L?^^-(4p+i-<i) _ (15b) 

The normalization of M(r) in Equ. (14) is chosen such that M = 1 for the 
extremal Reissner- Nordstrom solution. The limit a — > (in the original, unsealed 
variables) yields the generalised Schwarzschild solutions with constant M(r) = 
M, and thus 

2.1 New Variables — General 

The metric Ansatz (6) in terms of Schwarzschild coordinates used so far can 
not describe situations where the radius r of D — 1-spheres first increases and 
then decreases, as in the solutions studied in [27, 28, 29] with spaces of spherical 
topology. Even for asymptotically flat spaces with monotonically increasing r, 
the equations of motion resulting from the action (8) are singular when /i(r) = 0, 
i.e., at a horizon as well as for certain 'critical' solutions with a double zero of 
//(r). Following [2], we avoid this coordinate singularity by the most general 
static, spherically symmetric metric Ansatz 

ds^ = e^-(-)df - e2"(^)(ir2 - r\T)dnla_,) , (17) 

with r now a function of r, and substituting 

Adr^e''+^dr, ^ = {e-^ff , (18) 

with the notation x — dx/dr. 
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In order to obtain a system of first order differential equations (dynamical 
system), we introduce additional variables n, k, u, and Vp 



e r — n, 
re~^u — K — n, 



w — u 



re-^Hp^Vp , 



(19a) 
(19b) 
(19c) 
(19d) 



following the procedure in, e.g., [2, 30]. 

First we observe that the metric Ansatz (17) is explicitly invariant under 
reparametrisation of the new 'radial' variable r. Consequently varying the action 
cannot result in a differential equation for A, and thus all derivatives of e~^r in 
the action can be absorbed into a surface term that can be discarded. Then we 
introduce the new variables n, k, u, and Vp as Lagrange multipliers such that 
the resulting action contains at most one r-derivative and variation w.r.t. to the 
new variables yields Eqs. (19) as field equations. Finally we discard yet another 
surface term to obtain a compact expression. 

Each of these three steps is straightforward for p = 1, but requires some new 
techniques for p > 1 (see Appendix A for the details). As a result Eqs. (9) are 
expressed as 

d-2p-l 



S, 



I 



dr 



+{re~^n) — {k — n){e~^f 
Sm- j dTe^^'r'-'%d-2p-l){^^ 



2p 



[1 - n') 



(20a) 



1 WP 



+ w 



'Hi) 



^2 

2p 



u{e-^w) + ^{re-^Hp 
2p 



(20b) 



Varying the Action (20) w.r.t. A yields the 'reparametrization constraint' Ci — 
with 



Ci — {k — n)n 



d-2p 



2p 

(d-2p- 



"(1 



1 



V2(l-n2))p-i 
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1 2rr2\ 



2p 



(21) 



Varying w.r.t. n, k, m, and Vp by construction yields Eqs. (19) as field equations, 
and we will use them in the remainig variations. Varying w.r.t. A for fixed z/ + A 
yields 



re ^h — {K — n)n 



(r2(l -n2)) 
varying w.r.t. r for fixed e^^"^ j-d-'^v yields 



T/2 



2p- 



(22a) 
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/4p+l-d X, 

re K—( n — K (ft — n) 

\ p / 

^/ — 3,., (i — 2p — 1 _ 9,\ 
- 2(P - 1)t3^((^^"M - ^ )) 

+ ? (22b) 

and varying w.r.t. w and ifp yields the remaining field equations 

2{p-l)w— - + K + {d-Ap-l)n)u , (22c) 

re-^Vp = 2p{d - 2p - l)w^Hp - («; + (d - 4p)n) Vp , (22d) 

Next, we fix the freedom to reparametrise the radial variable r by the gauge 
choice — r and define Vp — {w'^ — iy~^v. 

2.2 New Variables — Specific 

The equations of motion resulting from the action (20) are singular when = 1, 
due to the term (3) with p > 1, or when = 1, due to the term (4) with 
p > 1. In order to avoid these singularities, wc have to introduce new variables 
specifically adapted to the form of the action (20), whereas the procedure used 
to obtain this action is essentially the same as used for other gravitating matter 
systems in [2, 7, 30, 31, 15]. The conical fixed point observed in [15] is caused by 
the new variables required to remove the singularity at — 1 from the equations 
of motion for that particular model. 

Here, we define new variables it = u/t, h = h/t, v = v/t where t = {w"^ — l)/r, 
and introduce an additional, redundant variable y — t^/{l — n^), resulting in the 
system of differential equations with polynomial r.h.s. 

f = rn, OT s = —sn , where s = , (23a) 

!> = K — n , (23b) 

n^{K-n)n-2yP{l-n'^)(u'^ + ^) , (23c) 



2p> 

/Ap+l-d \. . ^ r ^ fd-2p-l 
K = ( n -K){K-n) + 2f [(1 - n^) ( -) 

1 - 

^{p-\){\-Kn)({d-2p-\){—^w'y?)^u^^—\\ , (23d) 

u A - d — Ap—1 \ 
y = I {^^wu H uj 
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-2fn [(d - 2p - 1)(1 + wV) + i-a' + 2^)]}!, , 



w = {w^-l)u, (23f) 
tL^{d-2p-l)(l + (w^ - l)h'^)w - (2pwu + K + {d-Ap- 2)nju , (23g) 

(23h) 
(23i) 



h — v — {2pwu — n)h , 

V — 2p(d — 2p— l)w^h — (2pwu + k, + (d — Ap — l)rijv , 



subject to the two constraints 



where the expressions 

Ci — {k — n)n — (1 — n^) 



C2 = 0, 



d-2p-l 



2p 



2p 



-yP[(d-2p-l){^+w'h')-u'-^) 



obey 



Ci = -2/n( 



C2 — —r^y nCi + AwuC2 , 



(24a) 
(24b) 



(25a) 

(25b) 

(26a) 
(26b) 



and thus the two constraints are preserved by, and therefore compatible with the 
differential Eqs. (23). 

Note that Eqs. (23b-i) are independent of r (or s) and thus remain regular 
for r — > as well as for r — > 00, i.e., for s ^ 0. 



2.3 Singular Points of the Differential Equations 

The differential Eqs. (23) determine the derivatives of the dependent variables 
Y — (r, Z) (or Y — {s — r~^, Z)) where Z — {n, k, w, u, h, v, y) w.r.t. r, 

Y = /(y) . (27) 

They are singular when one of the dependent variables diverges and at the fixed 
points (f.p.s) of the Dynamical System (23). For each such f.p. Yq with /(lo) = 
we can introduce new variables Y — Y — Yq and linearise the equations 

Y = MY + 0(y2) . (28) 
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Excluding the possibility that the matrix M has eigenvalues with vanishing real 
part, we can rewrite Eq. (28) in terms of a suitable basis Y — {Y_,Y^) as 

with positive definite matrices M_ and M+. Due to the theory of dynamical 
systems there exists a 'stable manifold' Y = (Y_,F+ = 0{Y^)) of initial data 
such that y — > as r — > +00 as well as an 'unstable manifold' Y — (y_ = 
such that y — > as r — > —00. Prop. 1 of [7] states conditions 
such that these solutions can be characterised by and depend analytically on 
parameters determined at the f.p.; this however is not possible in general. 

First, there is the case u — > 00, occuring when w ±1 while w = ru remains 
finite. Such solutions are of no interest since Eqs. (19c) and (22c) exclude maxima 
of 1^1 with 1^1 > 1. 

Another type of singularity with n — > while ku ^ occurs at regular 
horizons and will be discussed below. 



2.3.1 Regular Origin 

A regular origin with r as r — —00 is described by the f.p. with n = k = 
w = 1. Eq. (23g) then yields either u = —{d — 2p — l)/2p or m = 1, but only the 
second of these solutions is compatible with a regular origin. 



2.3.1.1 RegulEir Origin without Higgs Field 

With /i = and 2p + 1 < d < 4p + 1 Eq. (23e) imphes 

yP^l, (30) 

and suggests the Ansatz 

w{r) = l-br^ + 0{r^) , (31a) 

u{r) = -2br + 0{r^) , (31b) 

n{r) = l- Cnr^ + 0{r^) , (31c) 

«;(r) = l + c«r2 + 0(r^) , (31d) 

where & is a free parameter, whereas 

Cn = c, = 26^ , (32) 

(compare [26]). Expressing Eqs. (23c-g) in terms of the dependent variables 
n — (1 — n)/r^, k — («; — l)/r^, w — (w — l)/r^, u — {u—l)/r, and y — {y — l)/r 
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as functions of r yields the equations 



^h- k + rffi , (33a) 

r— = {d-2){fi-k)+rh, (33b) 

r^^rfu,, (33c) 

r^ = -du + rfu , (33d) 
dr 

r^^-dy + rfy, (33e) 
or, using the hnear combinations n+ — {d — 2)n — k and n- — n — k, 

r^^r{{d-2)U-U), (33f) 

r^^-{d-3)n_ + r(U-U) . (33g) 

The nonhnear terms fn etc. are analytic functions of r and n etc., and thus 
Eqs. (33c-g) fulfill the assumptions of Prop. 1 of [7]. Furthermore Eqs. (25) can 
be written as 

^^k-n + rfc, ^^2n-Aw'' + rfc,, (34) 

and allow us to obtain the relation n(0) = k{0) = 2w^(0). This guarantees 
the local existence of a one parameter family of solutions with the boundary 
Conditions (31, 32), analytic in r and b, defined for all b and |r| < ^(6) with some 

m > 0. 

2.3.1.2 Regulcir Origin with Higgs Field 

We supplement Eqs. (31) with the Ansatz for the Higgs field, 

h(r) =ar + 0(r^) , i.e., h(r) = + 0(r^) , (35) 

2b 

such that $(a;, t) is regular near x — 0. In view of Eqs (23h,i) this requires d — Ap 
and consequently 

v(r) = (2p-l)ar + 0(r^) , i.e., v(r) = -(2p - 1)^ + 0{r^) . (36) 

2b 

Evaluating Eq. (23e) at the f.p. yields 

yP^{l + {2p-l)¥y\ (37) 
12 



We now have two free parameters, a and b, whereas Eq. (32) is replaced by 
c„ = 26^(1 + 7)Vf, c„ = C(7)c„, with C(7) = l-^^^ (38) 

and 7 = [2p — l)a^/46^. Expressing Eqs. (23c-i) in terms of the dependent 
variables n, k, w, and u as above together with = [2ph + ■u)/(4p — 1), /i_ = 
{{2p - l)h - v)/{{Ap - l)r), and y = {yP{l + 7) - l)/r with 7 = (2p - l)/^^ as 
functions of r yields the equations 



(4p-2)(C(7)n-/i)+r/«, (39b) 



^dr ^ ^^^^^ -i^ + T^fn, (39a) 

^ ^■^'^ ' ^^^^^ 

r — = -4p u + rfu , (39d) 
ar 



(ir 

dL 

dr 
dt 

r 



-iph.+rf~^_, (39f) 
-4py + rfy, (39g) 



or, rewriting Eqs (39a-b) in terms of the new variables n+ = (4p — 2)n — k and 

n_ = C(7)n - 

r^=r((4p-2)/^-/,), (39h) 

Furthermore Eqs. (25) can be written as 

^ = « - + r/c, , ^ = 2n(l + 7)"^^ - + r/c, . (40) 

We use the Constraints (24) to eliminate fi_ as well as Eq. (39i) with the /in- 
dependent 'eigenvalue'. In addition we obtain two relations between n(0), k(0), 
u'(O), and /i+(0). Eqs. (39c-h) then satisfy the assumptions of Prop. 1 of [7]. 
This guarantees the local existence of a two parameter family of solutions with 
the boundary Conditions (31, 35, 36, 38), analytic in r, a, and 6, and defined for 
all a, 6 7^ 0, and |r| < ^(a, h) with some ^(a, h) > 0. 
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2.3.2 Regulcir Horizons 

For regular, i.e., non- degenerate horizons we may use z — 1/k as the independent 
variable tending to zero (compare [31]), and thus replace Eqs. (23) by 

z-^T^--^z{l + zfr), (41a) 
az n 

zj-y = z{l + zU)y, (41b) 

z^r^z{l + zfr)r , (41c) 

z-^w = z{l + zfr)w , (41d) 
dz 

d _ _ 

z—u^-u + zfu, (41e) 

z-^h^z{l + zfr)ji , (41f) 
az 

z^v^-v + zfy, (41g) 

with expressions fr etc. determined from Eqs. (23) that are regular at 2; = 0, 
while n = 0{z) can be computed from the Constraint (24a). These equations 
satisfy the assumptions of Prop. 1 of [7] and thus guarantee the existence of a 
family of solutions with u ^ and v ^ as z ^ 0, and with finite limits t/j, 
Uh, Th, Wh, and hh for the remainig variables r, y, r, w, and h. Finally, the 
Constraint (24b) yields the relation yh = {w^ — lY /r\. 

In terms of the coordinate r the behaviour near the horizon is (performing a 
shift in r such that Th = 0) 

r(r)=r,(^l + ^j+0(r^), (42a) 

n{T)=nir + 0{t^) , (42b) 
k{t)=t'^ + 0{t) , (42c) 

w{r)^w, + ^-^^ + 0{r'), (42d) 

u{t)=WiT + 0{t^) , (42e) 

Mr)=ftJl-(p- 1)15^1 +^ + 0(r*), (42f) 



v(r)^hir + 0{r^) , (42g) 

with expressions rii, wi, and hi determined by the free parameters rh, Wh, and 
hh- 

We, thus, have a three parameter family of solutions that satisfy the black 
hole boundary Conditions (42) and are (except for the pole in k) analytic in 
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T, r/i, Wh, and hh, defined for all rh, Wh, hh, and |t| < ^{rh,Wh,hh) with some 
^{rh,Wh,hh) > 0. 

Solutions with a regular origin can be seen as the limits Vh ^ of those with 
a regular horizon. Consider a sequence of black hole solutions with finite limits 
hh/rh CLh and (1 — Wh)/r\ — > hh as Vh 0. Taken as functions of r they 
converge to a solution with regular origin and parameters a and h. Starting from 
the differential equations for w;(r) and h{r), derived form the action (9b), we 
can linearise around w = 1 and /i = 0, use a rescaled radial variable p = r/r/j, 
and neglect all terms that vanish as rh 0, resulting in two hypergeometric 
equations. We thus obtain the relations 

/ {p-l)d+p+l \ / V-3p+l \ 

bP^l ^-2p-i m, and abP-^^i ^p^^ UhK'' , (43) 

valid for d — Ap or, without Higgs field, for 2p + 1 < d < 4p + 1. Defining 

the ratios = b/bh and (a = O'/O'h) and evaluating the binomial coefficients in 
Eqs. (43) for p < 4 and d = Ap yields the values shown in Tab. 1, with the well 
known ratios Cb = ^ and Ca = 2 for p = 1 (compare Fig. 8b in [1]). 



Table 1: Ratios between boundary conditions for the limit r/^ — > of black 
monopoles and the corresponding regular monopoles for p = 1, 2, 3, and 4 in 
d — Ap spacetime dimensions. 



p 


(b 




1 


1 


2 


2 


1.91485422 


2.55313895 


3 


2.31107711 


2.81289292 


4 


2.60713329 


2.97958090 



2.3.3 Reissner-Nordstrom Fixed Point 

The Reissner-Nordstrom (RN) f.p. is characterised hyn = w = u = v = 0, 
h — a, r — 1, and k — kq where Kq — d — 2p — 1, as for the degenerate horizon 
of an extremal RN black hole. Introducing f — r — 1, u — ru, h — h — a, and 
K — K — Kq, we obtain 

f = n + fn , (44a) 
w^u , (44b) 
u = Kl{a'^ — 1) W — Kq U + fu , (44c) 

ti = v + fh, (44d) 



15 



V = -kq v + fy , (44e) 
h^Kon + f„, (44f) 
k^-2KoR + , (44g) 

with expressions fi = 0{{f,w,u,h,v,n,k.)^). Since f can be eliminated by the 
Constraint (24a), and assuming kq > 0, the hnearised Eqs. (44) have one unstable 
mode n with eigenvalue Kq, one stable mode k with eigenvalue —2k.q. The {h,v) 
subsystem contributes the stable mode v with eigenvalue — and the zero mode 
h + v/kq (that has to be considered as unstable due to the requirement /i — >■ 0). 
The {w,u) subsystem contributes two modes with eigenvalues 

A = -^(l± V4q;2-3) . (45) 

For oP" < 3/4 there are two stable oscillating modes with complex conjugate 
eigenvalues with negative real part, for 3/4 < < 1 there are two stable modes 
with negative eigenvalues, whereas for > 1 there is one stable and one unstable 
mode, i.e., one positive and one negative eigenvalue. 

Thus, the dimension of the stable manifold is 4 for < 1 and 3 for > 1, 
corresponding to a 3- resp. 2-parameter family of solutions converging to the f.p. 
as T ^ +00. However, it is possible to extend the 3-dimensional stable manifold 
for > 1 into the region 3/4 < < 1 and define a 3-dimensional submanifold 
of the 4-dimensional stable manifold by the requirement that the variables in 
Eqs. (44) decrease faster than e~^^'^l'^ . When a family of asymptotically flat 
solution reaches a critical limit with a double zero of /i(r), the interior part 
with r < 1 will be a member of the corresponding 2-parameter family, while the 
exterior part with r > 1 will be the exterior of the extremal RN solution with 

w = and h = a > ^J^J^. 

Note, however, that these solutions do not approach a degenerate horizon as 
r — cxo and r — » 1, because and hence A = /n diverge. They describe 
geodesically complete spacetimes, asymptotically like AdS2 x 5'°'"^, 

ds'' e^^°^dt^ - dr^ - dnla-2) ■ (46) 

2.3.4 Asymptotically Flat Infinity without Higgs Field 

Asymptotically fiat infinity with s ^ as r ^ H-oo is described by a f.p. with 
n = K = 1. In the absence of a Higgs field this requires w = ±1 and Eq. (23g) 
again yields either u = — (rf — 2p — l)/2p or m = 1, but this time only the first of 
these solutions is compatible with asymptotically flat infinity. 
Together with Eq. (14), this suggests the Ansatz 

771 

n{s) = l--z^ + 0{z^) , (47a) 
k{s)^1 + c^z^ + 0{z'') , (47b) 
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w{s)^l-cz^ + 0{z'^) , (47c) 
P 

with z — s(<^-2p-i)/2p^ where the total mass M — ^^±1=^777,? and c are free param- 
eters, while 

c. = ^l^^EzDm . (48) 

Expressing Eqs. (23c-g) in terms of the dependent variables n = (1 — n)/z'^, 
k— {hi — 1)/^^, w = {w — l)/z^, u = {u + {d — 2p — l)/2p)/z, and y = y/{szY 
as functions of s yields the equations 

3p + 1 — d 2 c 

s— = n + N, + z Jn , (49a) 

ds p 

dkd-Sp-l^^ Or /.r.! N 

s— = ^- n -K + z^ffi , 49b 

ds p 

s^^zU, (49c) 

, (49e) 



ds 



or, using the linear combinations n+ — n + k and fi_ = n{3p + 1 — d)/p + k, 

= (/. + /.), (49f) 

s-^ = n_ + z^[-^ U + fn) ■ 49g 

ds p ^ P ' 

The nonlinear terms /fi etc. are analytic functions of s, 2;, and n etc., and thus, 
Eqs. (49c-g) fulfill the assumptions of Prop. 1 of [7] with z as independent vari- 
able. Furthermore Eqs. (25) can be written as 

^ = + ^ = 2fiy-4*2 + zVc, , (50) 

z^ z^ 

and allow us to obtain the relations fi:(0) = n(0)((i — 3p — l)/p and n(0)y(0) = 
2iZ;^(0). This guarantees the local existence of a two parameter family of solutions 
with boundary Conditions (47), analytic in M, c, and ^, defined for all M, c, and 
|s| < ^(M, c) with some ^(M, c) > 0. 

2.3.5 Asymptotically Flat Infinity with Higgs Field 

Asymptotically flat infinity with Higgs field is described by the f.p. with n — k — 
1, w — u — V — and h — a > The difficulty here is, that w and u decrease 
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exponentially, approximately as Q-^/^^-^ar ^ while n, k, h, and v converge much 
more slowly, as some power of s = 1/r. 

Assuming solutions with finite energy, we can use the Definition (14) and the 
Constraints (25) to obtain expressions for n = r{n — 1) and k = r{n — 1), that 
are bounded as long as w, m, v, and the mass function M(r) are bounded. 
This suggests to derive differential equations for w, u, h — a, v, and M with r as 
independent variable 

^ = « + (51a) 
^ = (2p - l)a'w + fu , (51b) 

with nonlinear terms and /„ built from w, u, h — a, v, n, k, and s, while 
d{h — a)/dr, dv/dr, and 

^ - - 1)'"-' + - + «^ + $ . (510 

consist entirely of such nonlinear terms. This yields one dimensional stable and 
unstable manifolds together with a 'center' manifold. In order to further analyze 
the orbits in the center manifold, in this particular case simply w — u — 0, we 
intoduce h = r{h — a + v), v = rv and obtain 

s-—^sfM, (52a) 
ds 

'f,--h + 'fh, (52b) 

s^ = s/s, (52c) 

with nonlinear terms /m, fh, and fy that are bounded expressions in terms of 
M, h, and v. Consequently there exists a three parameter family of solutions, 
partially characterised by the total mass and the asymptotic value of v. 



3 Numerical Results 

3.1 Numerical Procedure 

We are mainly interested in regular monopole solutions, i.e., solutions of Eqs. (23) 
connecting a regular origin, r = 0, with asymptotically fiat infinity, r oo. In 
addition there are 'black monopolcs" starting from a non-degenerate horizon at 
r — Th- To better understand some limiting cases of these two types of solutions, 
we also need solutions ending at the RN fixed point w = 0, h = a a.t r— 1 (joined 
with the exterior extremal RN solution i(; = 0,/i = Q;forr>l). 
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We use a procedure that could be called 'shooting and matching'. For each 

of the four cases mentioned, we use a suitable independent variable < ^ < cxd. 
First we integrate Eqs. (23) from ^ = to some C, = ■Cmax with a slightly modified 
Runge-Kutta algorithm, using the shooting parameters a and b or hh and Wh 
to replace expressions that become undetermined at the singular starting point 
^ = 0. Next we express the solution at ^ = ^max in terms of stable and unstable 
modes suitable for the chosen endpoint ^ — > cxo. A solution will only converge to 
that endpoint for data on the 'stable manifold' where (the values of) the unstable 
modes are functions of the stable modes given by the solution of a system of 
integral equations. Finally, the shooting parameters must be adjusted to satisfy 
the matching condition that the solution at ^ = ^max lies on the stable manifold. 

As usual for nonlinear systems, the procedure outlined above has a limited 
domain of convergence and requires good approximate initial data. Since we will 
study families of solutions depending smoothly on one or two parameters, such 
approximate initial data are given by varying these parameters in small steps. 

In order to emphasise the similarities as well as the differences between the 
members of the hierarchy of p-EBPS solutions, we present results (figures and 
numbers) for p = 1, 2, 3, and 4, repeating results for p = 1 from [1, 2] (for l3 = 0, 
i.e., without Higgs potential). In addition we will use results obtained in [16] for 
the p-BPS hierarchy in flat space and in [26] for the p-Bartnik-McKinnon (BK) 
hierarchy. 

3.2 Regular Monopoles 

As noted above, the limit a ^ can be taken in two different ways. The 
'unsealed' variables f — ar and h — h/a satisfy a-dependent equations together 
with the a- independent boundary condition ^ — >■ 1 as f — > oo. For a — this 
yields the equations for the BPS monopole in flat space. Near f = we have 
w = 1 — br^ + 0{r'^) and h = df + O(f^) with d = b/ ^/2p — 1 dictated by the 
Bogomol'nyi equations, and b = 1/6, 0.1623, 0.13740, and 0.11775 for p = 1, 2, 
3, and 4 respectively (as found in [16] but normalised differntly). Starting from 
these values we could increase a and adjust the parameters d{a) and b{a) such 
that the boundary conditions h ^ 1 and -u; ^ as f — > cxo remain satisfied. 

Alternatively, we can use the 'rescaled' variables r and h with the a-dependent 
boundary condition /i — > a as r — > oo. Since there is no Higgs potential, the 
resulting equations are a-independent. For a — and thus h = they describe 
the EYM system [26]. Near r = we have w = 1 — br"^ + O(r^) with b = a% 
and h = ar + 0{r^) with a = a'^d, and thus a ^ b/ ^2p — 1 when a ^ 1. In the 
absence of a Higgs potential it suffices to adjust one of the parameters a or 6 such 
that as r —> oo, whereas a is determined by the solution as lim^^oo h{r). 

In addition to these fundamental monopole solutions with w{r) > 0, there 
exist excited solutions with N = 1, 2, . . . zeros of w{r). Starting from the N*^ 
generalised BK solutions with a — and b — b^^ as given in Tab. 2, we can 
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increase a and adjust b{a) such that u'(r) —> as r —> oo. 

Table 2: Parameters of generahsed Bartnik-McKinnon solutions for p = 1, 2, 3, 
and A in d — Ap spacetime dimensions. 



p 






M 


b? 


l,ip) 


1 


0.453716 


0.651726 


0.697040 


0.704878 


0.706420 


2 


0.415609 


0.532402 


0.574678 


0.588191 


0.593799 


3 


0.410745 


0.502893 


0.539835 


0.554239 


0.562700 


4 


0.410127 


0.489756 


0.523142 


0.537311 


0.547598 



For N ^ oo the BK solutions converge to a non-trivial limiting solution for 
r < 1, and to the exterior of an extremal RN black hole for r > 1. Starting from 
these limiting solutions with a — and b — b^ we obtain limiting monopole 
solutions. 

The properties of the excited and limiting solutions are dominated by the 
behaviour near the RN fixed point, where w decreases exponentially (in r) and 
oscillates with a frequency proportional to y/o^ 



solutions coverge to the limiting solution as 
that value. 



0.75. Consequently all excited 
0.75 and cease to exist beyond 





Figure 1: Fundamental solutions for p = 4 with a = 1.0, 1.6, cimax, 1-6, 1.3, 1.2, 
and 1.125; the dotted curves are for the monopole in flat space. 



We thus obtain smooth one parameter families of solutions. Some fundamen- 
tal solutions for p = 4 are shown in Fig. 1 and N — 2 excited ones for p = 2 are 
shown in Fig. 2. These families start with the (rescaled) flat monopole for = 
or the BK solutions for > and end with critical solutions where fi{r) has a 
double zero at r = 1 . These critical solutions consist of two parts describing two 
geodesically complete {t =const) spaces: a non-trivial interior part for r < 1 and 
a trivial part with w = 0, h = a, and A =const for r > 1, i.e., the exterior of an 
extremal RN black hole. The extremal RN black hole has a degenerate horizon 
at r = 1 and is a geodesically incomplete spacetime. The non-trivial interior part 
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Figure 2: Second excited solutions for p = 2 with a = 0.1, 0.15, 0.2, 0.35, 0.5 
and 0.8; the dotted curves are for the corresponding BK solution. 

describes a geodesically complete spacetime with A — > oo as r ^ 1. The initial 
data for these families are quite similar for p = 1, 2, 3, and 4 (see Fig. 3). 




Figure 3: Initial data for the regular monopole solutions for p = 2 and 4. 



An interesting structure emerges when we include the values of a obtained 
from these solutions. For N — and 1 < p < 4, a first increases from to some 

Omax and subsequently decreases to some Simultaneously the mass M starts 
as a times the mass Mflat = \^2p — 1/p of the flat space monopole, increases to a 
maximum Mmax > 1 and subsequently decreases to Mcr = 1 (see Figs. 4 and 5). 
In Tab. 3 we have collected some relevant numbers. The difference ccmax — is 
very small for p = 1, but quite large for p > 1. 

Table 3: Maximal and critical values of a and masses of the fundamental mono- 
pole solutions for p = 1, 2, 3, and 4. 



p 








1 


1.40303 


1.00022 


1.38585 


2 


1.43629 


1.01212 


1.18516 


3 


1.63363 


1.02921 


1.14018 


4 


1.83513 


1.04432 


1.12034 
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This may be a good point to comment on the numerical procedure described 
above. The possible values of ^max are limited severely by the necessity to suppress 
the unstable YM mode growing exponentially with ^ = r, whereas h{r) converges 
only slowly to its asymptotic value h{oo) — a. For the p — 4 solution with ccmax 
as shown in Tab. 3, h{r) differs by more than 10% for the largest possible values 
of ^max- Requiring that the unstable Higgs mode, h in Eq.(52b) vanishes at ^max 
yields excellent results for the BPS monopole in flat space. But in the gravitating 
case considered here, fj^ has slowly decreasing contributions proportional to n or 
K, preventing to estimate h{oo) with an error much better than 10~^. The results 
obtained with 'shooting and matching' as outlined, however, are accurate up 
to rounding and discretization errors, typically around 10~^° and are with that 
accuracy independent of the choice of ^max- 

The (first three) excited solutions for p = 1, 2, and 3 show a monotonic 
increase of a from to = y/0.75, with the mass increasing from its BK 
value to 1. For p = 4, however, the values of a first increase to some ftmax with 
^max > 1; thcu dccrcasc to some «min with Mjniii < 1, and finally increase to 
acr with M = 1 (see Figs. 4 and 5). One might speculate that higher excited 
solutions, with N > 3, exhibit such maxima and minima also for p < 4. 

The maxima of a (and M) observed for the fundamental monopoles give rise 
to two 'branches' of solutions. Furthermore, for a = Omax there exists a 'zero 
mode', indicating a change in the number of instabilities against small time- 
dependent perturbations. Whereas the 'lower' branch from to a^iax inherits the 
stability of the fiat space BPS monopole, the 'upper' branch from ctmax to has 
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0.5 1 1.5 a 0.5 1 1.5 a 



Figure 5: Mass of the regular monopole solutions with p — 1, 2, 3, and 4. 

one such instability. 

For ot <^ 1 the excited monopoles with mass a; Mbk + ^Maat may be seen 
as 'superposition' of a BK solution for ar <^ 1 and a fiat space monopole for 
r ^ 1, and inherit the instabilities of the BK solution. The N^^ BK solution for 
p = 1 has 2N such instabilities (A'^ radial and N sphaleronic ones) and this might 
equally be true for p > 1. The maxima and minima of a (and M) observed for 
the excited p = 4 monopoles give rise to three branches of solutions, where the 
'middle' branch from Omax to amin has one additional instability. 

3.3 Black Monopoles 

As noted in Subsect. 2.3.2 the limits r/j — > of black monopole solutions yields 
regular monopoles. Conversely the (fundamental, excited, or limiting) regular 
monopoles can be taken as starting points for corresponding black monopoles 
with Th <^ 1. 

Using 'rescaled' variables one can, e.g., choose a value and adjust one of the 
parameters hh or Wh such that w — >• as r — >• oo, whereas a. is again determined 
by the solution as limr^oo^(^)- We thus obtain smooth two parameter families 
of solutions. 

The initial data arc quite similar for different values of p and arc shown in 
Fig. 6 for p = 3 with Vh = 0.5 and for p = 4 with Vh = 1.5. For hh 1, and 
thus a <^ 1 they start with w/j ~ 1 for the fundamental, N = solution (tiny 
Schwarzschild black hole inside a fiat space BPS monopole), or with as for 
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0.9 0.92 0.94 0.96 0.98 1 0.2 0.4 0.6 0.8 1 

Figure 6: Initial data for the black monopole solutions with rh — 0.5 for p — 3 
and rh — 1-5 for p = 4. 




Figure 7: Excited, N = 2 black monopole solutions for p = 4 and a = 0.4, with 
Th = 0.2, 0.4, 0.6, 0.8, and Wh = 0.5, 0.4, 0.3, 0.2, and 0.1. The dashed lines 
are for the corresponding regular monopole. The dashed-dotted lines show the 
r/i-dependence of hh and Wh- 

the corresponding EYM black hole. A?" = 1, 2, . . ., or — > oo if < 1 (compare 
Fig. 3 of [7]). 

For Th < 1 the families end with critical solutions as for the regular monopoles. 
The critical fundamental solutions occur at a value acr{rh), starting for Vh = 
from ctcr as given in Tab. 3 for the regular monopoles and reaching the value 
some rh close to but less than 1. For larger rh the fundamental solutions 
converge to the limiting one as — > 0.75 and cease to exist beyond that value, 
as do the excited solutions for all r/j < 1. 

The situation is quite different for Vh > 1- All solutions end with = 
where they bifurcate with a regular, i.e., non-extremal RN black hole. Black 
monopole solutions with various values of rh are shown ior p — A, N — 2, a — 0.4 
in Fig. 7. 

The relatively simple domain of existence of the families of black monopoles 
described above gets rather more complicated by the existence of maxima and 
minima (of a for fixed Vh or vice versa), giving rise to various branches. Fig. 8 
shows these domains as arh fs. a for p = 1, 2, 3, and 4. Regular RN black holes 
only exist for r/i > 1 (above the dashed diagonal). Fundamental and excited 
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0.5 1 1.5 a 0.5 1 1.5 a 



Figure 8: Domains of existence for (fundamental and excited) black monopole 
solutions with p = 1, 2, 3, and 4. 

black monopoles always exist for r/j < 1 (below the diagonal) and < 0.75 (left 
of the dashed vertical line), but extend beyond that region. Exhibiting maxima 
and minima (dotted lines) they either end as critical solutions (solid lines for 

> 0.75, including a short vertical line near rh — 1) or bifurcate with a regular 
RN black hole (solid lines for < 0.75). 

The fundamental solutions with r/^ <^ 1 always have two branches, a 'lower' 
one from a = to some Q;max(^/i) and an 'upper' one from a-^g^{rh) to acr{rh)- 
Likewise, the fundamental and excited solutions with a <^ 1 always have two 
branches, a lower one from = to some ruiax{oi) and an upper one from from 
rmax(a) to rbif(«). 

For p = 1 these upper branches of the fundamental solution cease to exist 
near Vh — 1. For p > 1 they exist for all Th and are connected. In these cases 
continuity requires the existence of yet another branch near Vh and a — -\/0.75, 
shown as detail in Fig. 8. 

For p = 1 the upper branch of the excited solutions ceases to exist for larger 
values of a. For p > 1 the situation is somewhat more complicated. Increasing 
a, there first appears a minimum of Vh near the bifurcation, giving rise to a third 
branch. For p — 2 and 3 the maximum and minimum eventually meet, i.e., the 
upper branch disappears and the lower and third branch are joined (compare 
Fig. 9). For p = A, however, all three branches continue as 0, with the 

'upper' branch of the black monopoles corresponding to the 'middle' one of the 
regular monopoles. This correspondence suggests that the third branch observed 
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0.05 0.1 0.15 0.2 ar^ 0.25 0.05 0.1 0.15 0.2 ar^ 0.25 



Figure 9: Initial data for the N — 2 excited black monopole solutions for p = 3 
with a = 0.04, 0.05, . . ., 0.25. 

for p = 2 and 3 is related to the steep (almost vertical) increase of b vs. a shown 
in Fig. 4. 

4 Summary and Discussion 

When it comes to choosing a model for a gravitating monopole in higher dimen- 
sions, the number of options proliferate. Our choices are guided by the twin 
criteria of having a Bogomol'nyi lower bound that can be saturated in the flat 
space limit of the Yang-Mills-Higgs (YMH) subsystem, and that aside from the 
gravitational coupling and the Higgs VEV there feature no other dimensionful 
constants in the model. This restriction is made to simplify both the analytic 
and the numerical analyses. The YMH systems employed, in d = 4p spacetime, 
are those resulting from dimensional descent from i?^^"^ x which in fiat space 
support selfdual solutions [16]. The gravitational systems are then chosen to be 
those members of the gravitational hierarchy (4) with p = q, such that the rela- 
tion between the dimensions of the gravitational term in the Lagrangian (5) has 
exactly the same relation to the dimensions of the YMH term, for all p. 

We have carried out a precise quantitative analysis of gravitating static mono- 
poles, both regular and black, in ci = 4p spacetime dimensions. The models we 
have employed are in a sense very special, but the properties of their solutions are 
generic. Our template has been the gravitating Georgi-Glashow model in d — A, 
in the BPS limit, i.e., in the absence of the Higgs self interaction potential. This 
was the part of the subject of study in [1, 2] in which the Higgs potential was ab- 
sent. We have emulated the results of [1, 2] exactly here, and have thus achieved 
our aim of showing that the generic properties of gravitating monopoles are fully 
understood at least within the restricted choice of models that we have exercised. 
As such, the present study is a general preliminary investigation into the nature 
of gravitating monopoles in higher dimensions, both regular and black. 

In Section 2, we have presented an exhausive analytic analysis of the residual 
one dimensional action of the static system subject to spherical symmetry in 
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the spacelike dimensions. Especially prominent is the analysis of the Reissner- 
Nordstrom fixed point, and due to our restricted choice of systems, no conical fixed 
points [15] feature here. (Some numerical analysis of less restricted models was 
carried out which indicated the existence of conical fixed points, but we do not 
report on these here.) A preliminary step in the analysis for the gravitating Yang- 
Mills p— hierarchy was carried out before proceeding to the case of YMH ind — 4p. 
This underpins the numerical results of [26]. That was followed in Section 3 
by the numerical analysis of the 4p— dimensional gravitating monopoles, where 
the qualitative features discovered in [1, 2] for the p = 1 case were reproduced 
quantitatively for p = 1,2,3 with high accuracy. It is clear that the regular 
pattern observed modulo 4p up to p = 3 will repeat for all p. 

Having achieved our aim of exhibiting the regularity of features modulo 4p 
of the simplest class of gravitating monopoles, it is perhaps in order to point to 
natural succeeding investigations. In this context we would exclude the otherwise 
obvious choice of employinging the usual (p — 1) Einstein-Hilbert system in mod- 
els in d > 5, since this would result in the additional complication of encountering 
conical fixed points. As for introducing the usual Yang-Mills system in d > 5, 
this is excluded on the grounds that already in the fiat space limit it would cause 
the mass/energy to diverge because of the half-pure- gauge asymptotic decay of 
the YM connection in all dimensions. The following options are open. 

a) Extend the static spherically symmetric Ansatz (7) to allow for the electric 
YM connection Ao = u{r) ^ Hjo-i- These would yield the gravitating versions 
of the dyons constructed in [16]. 

b) Construct the gravitating monopoles in odd spacetime dimensions. The 
main difference of the YMH models to be employed in this task, and the ones used 
in the present work, is that unlike the latter the former arise from the dimensional 
descent from x where is now even and hence N > 2. As such, the 
gravity decoupling limits of these solutions do not saturate the Bogomol'nyi lower 
bound [22] , but this is not important since in any case the gravitating monopoles 
are not selfdual. In this respect gravitating monopoles in odd spacetimes are 
similar to those in even spacetimes for N > 3, with odd N, so that the two tasks 
should be performed in parallel. Another difference between = 1 and A^ > 1 
models is that in the A^ > 1 models expressions of the Lagrangians for increasing 
p become progressively more cumbersome, in contrast with the A^ = 1 models for 
which the p— hierarchy is quite uniform. For this reason in the A^ > 1 case it is 
reasonable to restrict to p = 2 to get a glimpse of the qualitative features. Were 
one to restrict attention to the p = 2 YMH systems, then presumably the most 
aesthetic (if not necessary) choice for the gravitational system would be the p = 2 
member of the gravitational hierarchy. Restricting to p = 2 YMH systems, the 
pertinent examples are d = 6 + 1, d = 5 + 1 and ci = 4 + 1 spacetimes. The p = 2, 
d — 3 + 1 case is also of some interest, as it presents very different properties 
from the usual p — 1, d = 3 + 1 case, notably supporting mutually attracting 
like-charged monopoles [32]. 
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A The New Variables 



With the substitutions Eqs. (18) the Action (9) takes the form 
Sg^-- [ dre^+'-^lr'-'P-'il-e- 



(53a) 



'M ■ 



I dr 6''+' r'-'^ W^-'{e-'^y + 



dw^2 d-2p-lWP 



dr 
dh 

2p^"' "dr 



2p 



1 / 



.xdHp^2 



+ (d-2p- l)w'H^ 



(53b) 



As a first step we want to absorb all r-derivatives of A in Eq. (53a) into a surface 
term. To write this in compact form we introduce a polynomial Fp uniquely 
determined by the property Fp{z) + 2zFp{z) = (1 — zy. We thus obtain 



Sg = J d^i^e-H)- j dre'^^ 



X ^d-2p-2 



(re-^z>)(e-V)Fp(e-'V) 



+ {d - 2p - 1) [—(1 - e-'V)P + {e-^f^Fpie 
Zp 



-2X.2. 



with 



$(x) = eV-'^'-^xFp(a;2) . 
Next, we consider an action of the form 



S = j dre^ L{(p,e~^(p) , 



(54) 
(55) 
(56) 



and introduce ipi as abbreviation for e ^(fi. Varying this action w.r.t. ipi yields 
the second order equations 



'£M„e-% = — -Y.g^^h, where M, 



^3 



(57) 



and we assume that the determinant \Mij\ does not vanish identically. Varying 
the action w.r.t. A yields the first order equation 



(58) 
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i.e., the reparametisation constraint. 

In order to obtain a system of first order differential equations we now want 
to introduce new variables ipi as Lagrange multipliers. For a Lagrangian of the 
form 

L(<^,^) = P{ip) + J2i^i^ljjKi,{cp) , (59) 

i.e., qudratic in the derivatives, this is achieved by the replacement 

L{ip,e'^Lp) ^ L = L -Y^{e~^ipi - ^i){e'^ipj - ipj) Kij 

= P + Y.ie-^iPii^j + Vie"Vj - . (60) 

We can generalise this procedure to the case with higher powers of derivatives by 
the replacement 

L{^, e-V) L{^, ^) + E(e-V^ - ■ (61) 

Varying the new action S = j' dr e^ L w.r.t. the new variables -0^ yields 

$:M,,.(e-V,--V',)=0, (62) 
j 

i.e., what used to be abbreviations are now field equations. Varying w.r.t. ipi 
or A and substituting these new field equations reproduces the previous results 
Eqs. (57) and (58). 

Applying this procedure to {Sq — J d^) in Eq. (54) and replacing $(e^'^f) by 
$(n) yields Eq. (20a). Applying the procedure to Eq. (53b) yields Eq. (20b). 
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